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, Abstract 

I , In brane-world cosmology gravitational waves can propagate in the 

■ higher dimensions (i.e., in the 'bulk'). In some appropriate regimes the 

bulk gravitational waves may be approximated by plane waves. We sys- 
£f) , tematically study five-dimensional gravitational waves that are algebra- 

^ ■ ically special and of type N. In the most physically relevant case the pro- 

' jected non-local stress tensor on the brane is formally equivalent to the 

energy-momentum tensor of a null fluid. Some exact solutions are studied 
to illustrate the features of these branes; in particular, we show explicity 
that any plane wave brane can be embedded into a 5-dimensional Sik- 
los spacetime. More importantly, it is possible that in some appropriate 
regime the bulk can be approximated by gravitational plane waves and 
thus may act as initial conditions for the gravitational field in the bulk 
(thereby enabling the field equations to be integrated on the brane) . 

a -1 

1 Introduction 

> ! „ 

• rH , In brane-world cosmology, matter fields and gauge interactions are confined 

' to a four-dimensional brane moving in a higher-dimensional "bulk" spacetime. 

This paradigm is motivated by string and M-theory; in particular, generalized 
Randall-Sundrum-type models pQ are relatively simple phenomenological five- 
dimensional (5D) models which capture some of the essential features of the 
dimensional reduction of Hof ava- Witten theory . In a recent analysis of the 
asymptotic dynamical evolution of perfect fluid brane-world cosmological models 
close to the initial singularity, it was found that for an appropriate range of the 
equation of state parameter an isotropic singularity is a past-attractor 0. It 
was subsequently argued that the initial cosmological singularity is isotropic in 
brane-world cosmological models. 
The 5D field equations are 

^G ab = -A 5 Wg ab + 4&T ab , A6 = -J. (1) 
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Here, a,b — 0, ...,4, A5 and g a b is the 5D cosmological constant and metric, 
respectively. The projected field equations on the brane are g|: 

2 k2 

G>„ = —Kg^v + k T^ v + 6— Sfiu — E^u , (2) 

where fx, v are brane indices (i.e. fx, v = 0, ..,3) A = A5/2 + A 2 k|/12 and 
k 2 = Ak|/6. The term is quadratic in and dominates at high energies 
(Too = P > A), and the five-dimensional Weyl tensor is felt on the brane via 
its projection, Z a \, — C ac bdn c n d , where n a is the unit normal vector to the 
brane. There may also be terms that arise from from 5D sources in the bulk 
other than the vacuum energy A5, such as a bulk dilaton field. In general, 
in the 4-dimensional picture the conservation equations do not determine all 
of the independent components of on the brane (and a complete higher- 
dimensional analysis, including the dynamics in the bulk, is necessary) [5]. 

In these models the gravitational field can also propagate in the extra di- 
mensions (i.e., in the 'bulk'). For example, there might occur thermal radiation 
of bulk gravitons [Hj- In particular, at sufficiently high energies particle inter- 
actions can produce 5D gravitons which are emitted into the bulk. Conversely, 
in models with a bulk black hole, there may be gravitational waves hitting the 
brane. At sufficiently large distances from the black hole these gravitational 
waves may be approximated as of type N. Alternatively, if the brane has low 
energy initially, energy can be transferred onto the brane by bulk particles such 
as gravitons; an equilibrium is expected to set in once the brane energy density 
reaches a limiting value. (For an alternative approach see, e.g., |7j.) In this 
paper we shall study the consequences of assuming that in some appropriate 
regimes the bulk gravitational waves can be approximated by plane waves. 



2 An example of a Type N bulk with a plane 
wave brane 

First, let us consider an example of an exact 5D solution of type N with a 
negative cosmological constant. This particular example will give us some hints 
what we can expect from the brane-world analysis. In particular, the example 
makes it possible to interpret some of the exact solutions on the brane discussed 
later. 

Consider the Siklos metric [B] (where we have dropped an overall scaling) 

ds 2 = (2dvdu + H(u, x, y, z)du 2 + dx 2 + dy 2 + dz 2 ) , z > 0. (3) 

This metric is an Einstein space (i.e., R a b — ^9ab) if the function H solves the 
equation: 

\dx 2 + dy 2 + dz 2 zdz) K ' ' y ' ' ( ) 

These solutions describe gravitational waves propagating in a negatively curved 
Einstein space. Hence, these metrics generalise the AdSs spaces to solutions of 
the Einstein equation with a negative cosmological constant containing gravita- 
tional waves. These are exactly the type of models we want to investigate for 
the bulk. 
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To see how these can generate anisotropic stresses in the bulk, we can cal- 
culate the Weyl tensor of the above solutions. In general the Weyl tensor has 
the following non- vanishing components (in a coordinate basis): 

1 d 2 
C uiuj - -^2Q^J H > 

1 / d 2 d 2 d 2 \ 
6z 2 V dx l dxi dx k J 

where (i,j,k) is a permutation of (x,y,z). Since £ a b — C a cbdn c n d , the Weyl 
tensor can induce anisotropic stresses on the brane. For example, the Weyl 
tensor corresponding to the gravitational mode H = x 2 — y 2 , yields the bulk 
stresses on the brane given by C uxux = —C uyuy = — -K. We also have oscillatory 
modes by choosing, for example, 

H(u, x, y, z) = (A + Cz i ){x 2 - y 2 ) cos(w + <p x ) + (D + Bz*) cos(u + ip 2 ), (6) 

where A, B, C, D, ipi and ip2 are arbitrary constants. We will later see that 
similar modes naturally enter the analysis of the brane. 

Let us, for the sake of illustration, consider a simple exact 5-dimensional 
solution of a brane of this type. We consider a case in which a brane is embedded 
in the 5-dimensional solution eq.JJJl. The brane is located at z — zq with normal 
unit vector n = zd/dz. Note that if 



dH 

dz 



0, (7) 



W6 obtain for th.6 6xtrinsic curv£Ltur6 oc Q^u-, where g^ v is the induced metric 
on the brane. Hence, we see from the junction conditions that there are some 
vacuum solutions 1 by suitably tuning the brane tension A in the usual way. The 
induced metric on the brane can be written 

ds brano = 2dwdu + H(u, x, y)du 2 + dx 2 + dy 2 , (8) 

and hence, is a plane wave. Using eq.QJ and the junction conditions, the func- 
tion H(u, x, y) = H(u, x, y, zq) obeys 



d 2 d 2 \ , , d 2 



(9) 



These brane solutions are therefore not the usual vacuum solutions; they are 
in general sourced by the projected Weyl tensor of the bulk (see e.g. [3]). The 
gravitational waves in the bulk propagate along the brane and induce, via the 
projected Weyl tensor, a non-local anisotropic stress. This stress mimics the 
effect of an electromagnetic wave travelling on the brane at the speed of light. 

We note that any plane wave brane can be embedded (locally) in a Siklos 
spacetime. Given H(u,x,y), a solution can be found by choosing 

d 2 d 2 

H(u,x iyi z) = D{z,V 2 2 )H(u,x,y), V 2 = — + — , (10) 



^^Also note that a propagating electromagnetic wave on the brane can support a non-zero 



OH 

dz 
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where D(^, V|) is the operator 

oo 

D(z,Vl) = X>(*)(Vl)\ (11) 

i=0 

and are defined iteratively: 

F {z) = l, F i+1 (z) = -£ ^(z)dz^dz. (12) 

The brane is located at zo and -ff is chosen such that eq.JTJ) is fulfilled. In the 
special case where H(u, x, y) is a polynomial in x and y, the above sum will 
terminate and H will only contain a finite number of terms. Furthermore, by 
choosing the function H(u,x,y) appropriately, we note that the above example 
can include branes of Bianchi type III, IV, V, VI^ and VII/j. The relationship 
between brane and bulk can be further studied in a full 5-dimensional setting. 



3 A general type N bulk 

Let us now be more systematic and assume that the 5D bulk is algebraically 
special and of type N. This puts a constraint on the 5D Weyl tensor which makes 
it possible to deduce the form of the non-local stresses from a brane point of 
view. For a 5D type N spacetime there exists a frame l a , fi a , m l a , i — 2,3,4 
such that dD]: 

C a bcd = iCiiiji^mltcm^. (13) 

The frame £ ai n a , m l a (the frame vector n a is not to be confused with the brane 
normal vector n°) is defined via the only non-zero contractions: 

£ a n a = £ a n a = 1, m i a m ia = 5 ij . 

The "electric part" of the Weyl tensor, £ a & = C ac bdn c n d , where n c is the normal 
vector on the brane, can for the type N bulk be written as 

£ ab = Cuu [L(™> c ) ~ mi(£ c n c )] [e b (min c ) - m{{l c n c )\ . (14) 

Furthermore, one can easily check that £ a bn b = £ a a = 0. Note that for a type 
N bulk we also have 

£abt b = 0. (15) 

This can be rewritten using the projection operator on the brane, g a b — g a b — 
n a n b , 

£ ac g c b £ b = £ ab l b = 0, £ b = g b c £ c . (16) 

Hence, the vector lb is the projection of the null vector £ b onto the brane. By 
contracting this vector with itself, we get 

l b i b = ~{tn a ) 2 . (17) 

The following analysis splits into two cases, according to whether £ a n a equals 
zero or not: 
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1. £ a n a = 0: l a = £ a and null. 



2. £ a n a ^ 0: £ a time-like. 

These cases have to be treated separately and have different interpretations 
on the brane. From a 5D point of view they are of the same type, but since 
the 5D spacetime is anisotropic the orientation of the brane with respect to £ a 
is of significance. More precisely, considering plane- wave spacetimes, the case 
£ a n a — corresponds to when the wave propagates parallel to the brane, and 
in the case l a n a ^ the wave hits the brane. 



3.1 The case £ a n a = 

Let us first investigate the consequence of £^ being a null vector. We note that 
E^v, the four-dimensional projected Weyl tensor on the brane, can be written 
in this case 

= - (J^) 4X, (is) 

where e is some appropriate function. Hence, this is formally equivalent to 
the energy-momentum tensor of a null fluid. Equivalently we can consider it 
as the energy-momentum tensor of an extreme tilted perfect fluid. Using a 
covariant decomposition of E^ v with respect to a preferred time-like vector 
being orthogonal to some 3-surface with metric , we have 



f — _ 
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U 



(19) 



The non-local energy terms are thus given by 

U = e(l u u v ) 2 , 
Q/j, = e{Lu v %, 

?V = el { J u) . (20) 

Here, angled brackets (■ • ■ ) denotes the projected, symmetric and tracefree part 
with respect to the metric of the spatial 3-surfaces. The equations on the 
brane now close and the dynamical behaviour can be analysed (see later, and 
[TT]). Note that 

UV liV = Q lx Q v -^g lxv Q x Q x , (21) 

so that in this case is determined completely by U and Q M . 



3.2 The case tn a ^ 

In this case the vector l a has a component orthogonal to the brane. This implies 
that the vector £ a is time-like. This vector lives on the brane, and hence we can 
set parallel to £ fJ -. In this frame, the requirement l^Z^v = implies that we 
have U = = Q M and hence, we can write 

<W = - (±) V„ v . (22) 
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This is formally equivalent to a fluid which possesses anisotropic stresses with no 
energy density or energy-flux. From a brane point of view it appears as if these 
stresses are superluminal; however, as can be seen from from a 5D point of view, 
this is just an artifact of living on a brane in a higher-dimensional spacetimc. 
The stresses do have a gravitational origin, namely from gravitational waves in 
the 5D bulk. 



3.3 Exact solutions describing radiating branes 

Next, we will find some exact solutions in the case £ a n a ^ 0, which are very 
illustrative and describe features of branes that seems to have gone unnoticed 
in the literature. We assume that the brane contains an isotropic fluid, i.e. the 
energy-momentum tensor takes the form 

T^u = {p + P)u^,u v + pg^, (23) 

and that p = (7 — l)p, = U = 0. We also assume that the brane has isotropic 
curvature, T^u v ) — (e.g., is of Bianchi type I or V). From the propagation 
equation for the non-local dark energy U, the simple constraint 

o-^V^ = 0, (24) 

where a^ LV = V \^u v ^ is the shear, is obtained (see 52)- Since both cr M „ and 
V^u are symmetric, space-like, trace-free and of maximal rank 3, and defining 
the natural inner product, the constraint u^ v V^ v = implies that u^ v and V^ v 
are orthogonal. We have the evolution equation 



where 9 = V''u M and £l a v is antisymmetric and is the rotation tensor with 
respect to a set of Fermi-propagated axes (recall that an overdot is defined by 
' = u^Vp). Writing 9 — 3a/ a, the shear equations then give 

d t {a^) + Q- a a^ = Q. (26) 

This can easily be integrated to give cr fJiV cr IJ ' v = A 2 a~ e . Thus the shear scalar 
decays withe the usual a~ e behaviour and is unaffected by the non-zero V^ v . 
This is simply due to the fact that V^v is orthogonal to a^ v and thus it cannot 
affect the shear scalar. The generalised Friedmann equation now becomes ^2] 

£! = A + + (■> + 1_p°_\ tor) 

a 2 3 3 a 6 6 a 2 3 \ 2A J ' [ ' 

where, for example, k = for Bianchi type I and k = — 1 for Bianchi type V. This 
equation can now easily be solved in quadrature. Due to the quadratic terms in 
the energy-densities arising from brane effects, the shear mode is subdominant 
in the past as long as the isotropic fluid is stiffer than dust. 
It is convenient to define two new variables 
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In tensor form, the shear evolution equations are 



<9 t X+ [n,X] =uN, 



(28) 



where ui = 6/(k 2 A). This equation can be integrated in quadrature once we 
know Y (recall that from a brane point of view, there are no evolution equations 
for Y). We solve the above equations by choosing the particular gauge in which 
f2 = (X is not diagonal in general) , so that the axes are Fermi-propagated. A 
solution can now be written: 



In this solution it is easy to show that if Xo, Yo are orthogonal initially, then X 
and Y will be orthogonal at all times. From the assumption of a periodic Y, it 
follows that the function H satisfies a wave equation as desired. 

The above solutions describe brane-worlds where the shear oscillates with a 
frequency 



This oscillation is driven by a similar oscillation in the non-local anisotropic 
stresses "P M „. The frequency of the oscillation is given by the brane-tension and 
the gravitational constant on the brane and hence the physics on the brane. 
If these gravitational waves came from a source in the bulk, then that source 
must somehow know of the physics on the far-away brane. Unless there is a 
very special configuration in the bulk, then this situation does not seem very 
likely. A more likely scenario is that it is the brane itself which is the source of 
the gravitational radiation. The brane oscillates with a given frequency which 
is dictated by the brane tension and the gravitational constant on the brane. 
This oscillation causes the brane to emit gravitational radiation into the bulk. 
The frequency of the gravitational radiation will thus carry a distinct signature 
given by its frequency. If the bulk consists solely of an incoming or an outgoing 
wave then the type N approximation is valid. Suppose there is a source in the 
bulk emitting radiation with a different characteristic frequency than the brane. 
Then close to the membrane there will be a mixture of incoming and outgoing 
waves, unless some sort of equilibrium is established. This implies that the 
simple assumption that there exists a preferred wave vector just outside the 
brane (and the bulk just outside the brane is of type N) breaks down. 

4 The effect of gravitational waves on the dy- 
namical behaviour of the brane 

Finally, we will investigate what kind of effect this type N bulk may have on the 
cosmological evolution of the brane. It is believed that an isotropic singularity 
is the most likely initial starting point for a classical brane- world [3]. For this to 
be correct, the initial singularity has to be stable into the past. Let us assume, 
therefore, that we are in the regime of an isotropic past and the cosmological 
evolution is dominated by an isotropic perfect fluid with equation of state p = 



X = X cos ujt + Y sin wt, 
Y = Yq cosujt — Xo sinwi. 



(29) 
(30) 
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(7 — In this approximation, and assuming the case £ a n a = 0, the equations 
for U and Q M are 

U + ~M = 0, ^ + -^ = 0. (31) 

For the isotropic singularity, the expansion factor is given by 9 = l/(^t). Hence, 
defining the expansion-normalised non-local density, U = U/9 2 , and non-local 
energy flux, = Q^/9 2 we get 

U=^fry-2)U, g Al -^(3 7 -2)Q M . (32) 

Hence, the isotropic singularity is stable to the past with regards to these stresses 
if 7 > 2/3. 

Similarly, for the case £ a n a = we have a FRW universe to the future with 
9 = 2/(7*), and thus 

^=|-(37-4)C/, M = T-(37-4)Q M . (33) 
3 7 3 7 

This implies that if the isotropic fluid on the brane is stiffer than radiation 
(7 = 4/3), then the FRW universe is unstable to the future with respect to 
the non-local stresses. However, since dust has 7=1, our physical universe is 
believed to be stable with respect to these non-local stresses. 



5 Discussion 

These phenomenological results are further supported by a more detailed anal- 
ysis of the asymptotic behaviour of two tilting 7-law fluids in a class of Bianchi 
type VIo models I n particular, the physically relevant case of interest here, 
namely that the second fluid is a null fluid or fluid with extreme tilt, is investi- 
gated. All equilibrium points are found and their stability determined, so that 
the local attractors can be established. It is found that the dynamical effects of 
the projected Weyl tensor are not significant asymptotically at early and late 
times. 

In brane cosmology the only degrees of freedom in the bulk are the higher- 
dimensional gravitational waves, which propagate in the bulk The prob- 
lem of initial conditions for these bulk gravitational waves are of great impor- 
tance |13j . Indeed, we need initial conditions for the gravitational field in the 
bulk to be able to integrate the field equations and determine the dynamics on 
the brane. We have argued that in some situations the bulk gravitational waves 
can be approximated as type N plane- wave solutions. Consequently, perhaps an 
appropriate set of initial conditions is to assume plane behaviour in a suitable 
regime of the bulk. This can be achieved by demanding that the 5-dimensional 
bulk is algebraically special and of type N as above. The effects this type N 
bulk may have to the cosmological evolution of the brane is studied in JT] . 
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